Molecular Genetics and Genomics

Codon usage bias is a fundamental genomic characteristic that prefers non-random preferential use of synonymous codons. It is a major determinant of translational efficiency, gene regulation, and molecular evolution. However, the evolutionary bias and functional relevance of codon usage bias across the plant lineage is poorly defined and yet to understand what are the major factors responsible for relative synonymous codon usage (RSCU) in genomes and how codon usage bias influences the gene regulation, molecular evolution genomes. A genome-wide codon usage bias study of coding DNA sequences of 262 plant genome was conducted. It encompassed more than 4.6 billion codons from > 11 million coding sequences. Relative synonymous codon usage, codon adaptation index, codon-anticodon mapping, effective number of codon (ENC)-GC3, GC1,2-GC3, parity rule 2 (PR2-bias), molecular economy, and machine learning approaches were used for the study. It was found that codon usage bias was strongly non-random and exhibited a clear phylogenetic structuring. The higher plants favoured A/T-ending, whereas early-diverging lineages were enriched in G/C-ending codons. Analysis of RSCU, codon adaptation index, and codon-anticodon pairing indicated that translational selection is mediated by tRNA availability, contributing sustainability to these molecular patterns. Machine-learning approaches identified a small subset of codons having outsized influence on genome-wide codon usage landscapes. Further studies revealed the presence of robust inverse relationships between the effective number of codons and GC content at synonymous third positions. Neutrality analysis revealed approximately 61% of variation was driven by mutational pressure, tempered by selective constraints. Phylogenetic reconstruction showed a progressive relaxation of codon bias from algae to angiosperms while maintaining a conserved molecular economy cost of ~ 30 ATP per codon across the lineages. The study revealed codon usage bias is lineage-specific evolutionary conserved trait governed by mutation, selection, and translational optimization.

RNA interference (RNAi) is a widely exploited reverse-genetics tool with potential uses for disease control. Successful RNAi has been reported in trematode-vectoring snails, but the composition of RNAi effector-encoding gene complements, a key driver for RNAi efficiency, remain unstudied in these species. Using bioinformatics and comparative genomics, we searched for orthologues of 115 RNAi effector sequences in genomes or transcriptomes of four snail vectors: Biomphalaria glabrata, B. pfeifferi, Bulinus truncatus, and Lymnaea staginalis. Gene expression patterns of selected RNAi effectors were then examined across developmental stages and tissues of the model B. glabrata snail. At least 74 RNAi-related proteins were conserved across all four species, including core components known to be essential for gene silencing. Classical systemic RNAi-deficient (SID) genes that facilitate systemic RNAi in other systems were absent, suggesting that alternative pathways may compensate for dsRNA uptake and transport. Core effectors of secondary RNAi amplification and heritable RNAi were not detected. Expressions of Dicer-1, Argonaute-2, and the exonuclease Eri-1 did not vary significantly with snail size. A putative RNAi-inhibiting Staufen orthologue showed elevated expression in the ovotestis, while another putative cholesterol-interacting gene was overexpressed in the trunk tissue and may partly contribute to RNAi import. Altogether, our results present the most comprehensive overview of RNAi pathway effectors in major intermediate snail hosts for trematodes. The findings underscore the likely broad potential for RNAi use in trematode intermediate hosts as an experimental tool and potential control method.

BackgroundThe corn smut fungus Ustilago maydis is an important microbial model organism representing a genetically amenable and readily cultivable basidiomycete. Research in this fungus addresses a broad range of fundamental questions and its biotechnological exploitation is on the rise. Although genetic engineering in principle is well established, efficient methodology for synthetic biology approaches such as metabolic engineering or pathway transplantation has remained limited. ResultsHere, we present a comprehensive toolbox for U. maydis based on modular cloning and the characterization of more than 20 promoters. Careful comparative evaluation of insertion loci and terminator as well as reporter effects was conducted and a novel color-based strategy for straightforward genome integration was implemented. Moreover, the cloning and subsequent one-step integration of four transcriptional units into U. maydis was demonstrated by creating a "rainbow" strain producing four fluorescent proteins. ConclusionOverall, this next generation toolkit strongly advances genetic engineering and systems biology approaches in U. maydis, fostering its development into a valuable and competitive fungal chassis and prime model, particularly in applied research.

Demographic factors are intrinsically crucial to evaluate species' extinction risk. However, measuring them remains difficult and time-consuming and the use of genomic summary statistics has been advocated to assess the conservation status of a species. In the present study, we estimated (i) the census number (Nc), (ii) effective population size (Ne) over three different time periods, recent, historical and ancient, (iii) neutral genetic diversity ({pi}4), and (iv) a measure of the efficacy of purifying selection ({pi}0/{pi}4) for 101 plant species using population genomic sequencing data. Twenty-one species are from the Plant Species with Extremely Small Populations (PSESP) program of SW China. Threatened species exhibited significantly lower Ne, Nc, {pi}4, and weaker purifying selection, but had a higher Ne/Nc ratio than non-threatened ones. Nc was the main determinant in identifying conservation status, and contemporary neutral genetic diversity was predominantly influenced by historical Ne. In the absence of demographic information, genetic parameters are a good proxy of conservation status, likely because currently threatened species also had a low historical population size. In summary, our findings suggest that direct estimates of Nc are more useful than {pi}4, although the latter remains a valuable conservation indicator. Hence, efforts such as the PSESP should be extended.

Developmental processes depend on carefully coordinated gene expression. Expression is modulated by the binding of transcription factors (TFs) to cis-regulatory elements (CREs), like enhancers and promoters. Many computational and experimental approaches have been developed to find CREs, particularly enhancers, in the genome, each with strengths and caveats. Given the increasing availability of ATAC-seq data and methods to find TF binding therein, we hypothesized that we could use TF footprinting tools to find clusters of TF binding events within accessible chromatin that may act as CREs. Using Drosophila anterior-posterior patterning network as a test bed, we used a digital genomic footprinting tool (DGT), TOBIAS, on previously published early embryo ATAC-seq data to characterize the TF footprint landscape of 16 TFs essential for embryonic patterning. Even in this system, with its extensive enhancer annotation, most footprinted TF binding sites lie outside of known enhancers, with intergenic and intronic regions hosting the highest TF footprint count, albeit at low density. To find potential novel enhancers, we identified high-density TF footprint clusters that are highly conserved and overlap with active enhancer histone mark signals. Five high confidence candidates were selected for reporter assay validation and all five were found to drive spatially patterned expression in the embryo. This study shows that even in a highly characterized system, the analysis of footprinted TF binding sites in ATAC-seq data can uncover new regulatory regions and suggests this approach may be helpful in using existing ATAC-seq data to find novel CREs. ARTICLE SUMMARYGiven the increasing availability of ATAC-seq datasets, workflows to exploit the data to uncover new cis-regulatory elements (CREs), including enhancers, are valuable. Using early anterior-posterior patterning in the Drosophila embryo as a test case, we find that previously published transcription factor footprinting tools and ATAC-seq data can be analyzed to yield new candidate CREs. Experimental validation confirms the activity of selected candidate CREs, suggesting that existing data can be analyzed to find novel regulatory elements.

Congenital disorders of glycosylation and deglycosylation are rare, serious, and lethal disorders afflicting humans. CDGs and CDDGs result in loss of function enzymes which fail to build or break down oligosaccharides on proteins. This can produce protein aggregates and, in turn, reactive oxygen species that harm the cell eventually leading to autophagy and apoptosis. Because sperm contain high concentrations of polyunsaturated fatty acids, they are especially sensitive to these effects, which is understood as one of the leading factors in human male infertility. Sperm are developed in zebrafish similarly to humans and are useful models to examine human reproductive health, as well as genetic disorders. The combination of these advantages makes the analysis of sperm from zebrafish with heterozygous ALG1 or DPAGT1 CDGs or the NGLY1 CDDG suitable. Analysis of sperm concentration, motility, status, viability, and hypoosmotic swelling demonstrated the effects of these disorders on sperm quality. Results showed a significant decrease in sperm concentration, motility, and hypoosmotic swelling for all mutant zebrafish compared to the wild type. This suggests that CDGs and CDDGs influence the amount of sperm produced, the percentage of sperm cells that are mobile, and the integrity of the plasma membrane.

Gene model for the ortholog of raptor in the D. grimshawi May 2011 (Agencourt dgri_caf1/DgriCAF1) Genome Assembly (GenBank Accession: GCA_000005155.1) of Drosophila grimshawi. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of raptor in the May 2011 (Agencourt Dere_CAF1/DereCAF1) Genome Assembly (GenBank Accession: GCA_000005135.1) of Drosophila erecta. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Density regulated protein (DENR) in the Apr. 2013 (BCM-HGSC/Deug_2.0) (DeugGB2) Genome Assembly (GenBank Accession: GCA_000236325.2) of D. eugracilis. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Ectomycorrhizal (EM) fungi are vital to forest ecosystems, supporting tree growth and survival. However, their inclusion in conservation policy and action remains limited and little is known about the status of their genetic diversity, which is essential for their long-term survival and adaptation. The Global Biodiversity Framework adopted a genetic indicator based on the effective population size, Ne, to monitor genetic diversity in all species. To date, it is still uncertain how Ne, a key parameter, can be reliably assessed in species with complex life history traits. Ectomycorrhizal fungi are a highly diverse group of taxa displaying haplodiplontic life cycles with partially clonal reproduction. Here, we review the literature to understand how these life history traits might affect Ne and its estimation in six species of EM fungi. We estimated Ne in 19 populations using eight genetic and genomic datasets from selected studies. We compared Ne estimates using Linkage Disequilibrium (LD) and Sibship Frequency (SF) methods. We tested how Ne estimates change due to partial clonality and genetic structure gradients and whether the number of genetic markers influence the precision of the estimates. We show a systematic bias in Ne estimations when large clones are present and when populations are not correctly delimited. We found both methods are not robust to these factors, which makes them unreliable for conservation assessment purposes in EM fungi. This study provides new perspectives for further research into the links between life history traits and the effective population size of ectomycorrhizal fungi.

Gene model for the ortholog of tango (tgo) in the Sep. 2015 (UC Berkeley ASM127793v1/DbusGB1) Genome Assembly (GenBank Accession: GCA_001277935.1) of Drosophila busckii. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Organisms have spent their life histories exposed to background levels of natural ionizing radiation. To document the role that radiation plays, the deprivation of these natural levels has been studied by incubating organisms in the shielded space of underground laboratories. We report here on two studies (Study I and Study II) using Aedes aegypti for the first time as a model organism incubated 655 meters underground at the Waste Isolation Pilot Plant (WIPP) outside of Carlsbad, New Mexico, U.S.A. Male mosquitos were incubated at the surface exposed to natural background radiation, and were compared to two underground treatments in which incubators were supplemented with radiation sources used to mimic background and these groups were compared to the underground, radiation-deprived treatment. In Study I, the mosquitos incubated underground in the absence of natural radiation had higher levels of mortality compared to those incubated at the surface and PCA plots of the two transcriptomes were clearly differentiated. Study II was conducted the following year and the experiment was narrowed to include only the surface control and underground, radiation-deprived treatment which allowed for four biological replicates. Again, there was a higher level of mortality in the mosquitos grown underground compared mosquitos grown at the surface. Transcriptomes were not as clearly differentiated by PCA analysis and fecundity data were similar between the two groups. Functional analysis of transcriptomic DEGs from two independent studies suggested there are stress responses in radiation deprived mosquitoes. The absence of a secondary stressor in Study II is discussed as an explanation for the transcriptome differences in the two experiments.

Anoxia is a major stress for most vertebrates and frequently accompanies harsh winter conditions, particularly in species that spend much of the season frozen solid. North American freeze-tolerant wood frogs (Rana sylvatica) can survive several months without oxygen and endure whole-body freezing for up to eight months of the year, with [~]70% of total body water frozen as extracellular ice, yet revive when temperatures rise in spring. Survival depends on multiple adaptations, including tolerance of prolonged oxygen deprivation while frozen, when breathing and circulation are halted. A key strategy involves hepatic glycogen mobilization, producing large amounts of glucose that are distributed to tissues where it functions both as a cryoprotectant and as a substrate for anaerobic ATP production. The present study examines the role of histone lysine methylation and demethylation in regulating liver proteins under anoxic conditions. Relative protein expression of seven histone methyltransferases (ASH2L-S, ASH2L-L, RBBP5, SETD8, SMYD2, ESET, SETD1), six lysine demethylases (KDM1A, KDM3B, KDM4A, KDM4B, KDM5A, KDM5C), and eight histone marks (H3K4me1, H3K4me2, H3K9me3, H3K27me3, H3K36me3, H3K79me3, H4K20me1, H4K20me3) were evaluated in wood frog liver under control, 4-hour, and 24-hour anoxia exposures. The data indicate that histone lysine methylation and demethylation contribute significantly to transcriptional regulation under anoxia. Specifically, H3K4, H3K36, and H3K79 methylation were associated with transcriptional activation, whereas H3K9, H3K27, and H4K20 methylation correlated with transcriptional repression. These findings highlight the dynamic role of epigenetic regulation in supporting hypometabolism and stress adaptation in freeze-tolerant wood frogs.

BackgroundHeterochromatic long arm of mouse Y chromosome harbors the multicopy species-specific sequences Ssty, Sly, Asty and Orly that are transcribed in testis and have known functions in male fertility. Of these Ssty and Sly encode proteins - yet all the transcripts are not translated. To investigate the roles of these Y-heterochromatic transcripts further, we analyzed them. MethodsMice with 2/3rd deletion of the Y-chromosome (XYRIIIqdel) and its wild type (XYRIII) were used in this study. Bioinformatic approaches, small RNA northern blots, Electrophoretic Mobility Shift Assays, Luciferase reporter assays, dPCR analysis, RT-qPCR assays and western blotting techniques were used to identify piRNAs that regulate autosomal genes. ResultsWe demonstrate that the multicopy gene families from mouse Y-long arm generate piRNAs predominantly in testis. We observed sequences homologous to these piRNAs in the UTRs of a few autosomal genes, which are differentially expressed in the sperms of XYRIIIqdel mice. Furthermore, the Endogenous Retrovirus Element (ERV) LTR, found in the Orly1 transcript identified piRNAs in the database, showed homology to UTRs and associated genomic regions of a few autosomal genes. Orly1 showed a reduction in genomic copy number by digital PCR in XYRIIIqdel mice. One of the four autosomal genes containing the ERV segment in their UTRs, showed a differential testicular protein expression in the mutant mice. ConclusionsThus, we further elucidate that different classes of repeats from Y-chromosome regulate autosomal gene expression via piRNAs. Besides, this study also identified novel roles for a Y-derived ERV in autosomal gene regulation in testis.

ObjectiveTo develop and validate a cell-free DNA (cfDNA) fragmentomic classifier for the early prediction of spontaneous preterm birth (PTB) using routine first-trimester non-invasive prenatal testing (NIPT) data. MethodsA nested case-control study was conducted within a prospective multicenter Vietnamese cohort comprising 286 pregnancies, including 82 spontaneous PTB cases and 204 term controls. Maternal plasma cfDNA collected during routine first-trimester NIPT (median gestational age, 12 weeks) was sequenced to a depth of approximately 20 million reads per sample. Five fragmentomic feature categories including copy number alterations, end-motif composition, nucleosome distance, fragment length, and joint fragment-lengthxend-motif were evaluated for PTB prediction. Machine learning classifiers were developed in a training cohort (n = 228, 65 PTB vs 163TB) and tested in a validation cohort (n = 58, 17 PTB vs 41 TB). ResultsAmong the five fragmentomic feature classes evaluated, 4-mer end-motif (EM) profiles exhibited the most pronounced differences between PTB and term control samples. Consistent with these findings, the EM-based classifier demonstrated the highest discriminative performance in the validation cohort, achieving an AUC of 0.970 (95% CI, 0.912-1.000). At a specificity >90%, the model achieved a sensitivity of 94% (95% CI, 78-100%). ConclusionThese findings demonstrate that cfDNA EM signatures derived from routine first-trimester NIPT can accurately identify pregnancies at risk of spontaneous preterm birth, without additional blood collection or sequencing, thereby extending the clinical utility of existing prenatal screening infrastructure. KEY POINTSO_ST_ABSWhat is already known about this topic?C_ST_ABSO_LICurrent first-trimester prediction strategies based on maternal characteristics, cervical length, and biochemical markers have limited predictive accuracy, particularly in nulliparous women. C_LIO_LIExisting cfDNA-based approaches have shown only modest performance or require additional assays, limiting clinical applicability. C_LI What does this study add?O_LIExisting NIPT sequencing data can be repurposed (without additional blood sampling or sequencing) for accurate prediction of spontaneous preterm birth (AUC=0.970). C_LIO_LIA classifier employing 4-mer end-motif (EM) profiles achieved an AUC of 0.970. At a specificity >90%, the model achieved a sensitivity of 94%. C_LI

The NAD+-dependent histone deacetylase Sirt6 regulates transcription of multiple classes of genes, including those involved in metabolism, immune response, oxidative stress response, and development. Defining the Sirt6-regulated transcriptome is relevant to understanding the various important physiological roles of Sirt6, such as extending lifespan, maintaining metabolic health, and tumor suppression. Numerous studies have identified Sirt6 target genes, using both targeted and genome-wide approaches; however, consensus is limited and there has yet to be a systematic analysis of gene expression changes induced by altering Sirt6 levels. In the present study, we conducted a meta-analysis of 19 mammalian RNA-seq datasets in which Sirt6 levels were perturbed (knockout, knockdown, or overexpression). Our analyses included Gene Set Enrichment Analysis, pathway analysis of differentially expressed genes, and identification of individual differentially expressed genes. Our analysis identified consistent gene expression changes associated with lowering Sirt6 levels, including increased expression of immune response and ribosomal protein genes and reduced expression of lipid oxidation and oxidative phosphorylation genes. Extracellular Matrix and E2F target genes also had consistently increased expression upon Sirt6 reduction, highlighting novel regulation by Sirt6. To determine the conservation of gene regulation by Sirt6, we performed additional RNA-Seq meta-analysis on tissues from Drosophila melanogaster with Sirt6 deletion and overexpression. The fly datasets produced similar results to the mammal results, except for lipid oxidation genes, which were found to increase in Sirt6-low conditions. These results provide consensus about conserved and novel pathways transcriptionally regulated by Sirt6.

Common variants detected by the genome-wide association studies (GWAS) create a wealth of knowledge on genetic component of individual traits and diseases. Elucidating the molecular mechanism behind the vast majority of these variants that are found to be non-coding remains a largely unsolved task, especially when distal and pleiotropic interactions between regulatory elements where these variants occur and gene promoters are taken into account. Focusing on four diseases with immune-mediated mechanisms namely ulcerative colitis, Crohn's disease, primary sclerosing cholangitis and ankylosing spondylitis, we demonstrate the utility of the targPred tool, providing prediction of genes targeted by the regulatory variants. We demonstrate that taking into account evolutionary and comparative genomic data, previously unobserved mechanistic trends (the platelet, vascular and sterol clusters) can be detected in terms of implicated genes targeted by the regulatory elements containing common variants, shared between all four diseases, as well as specific trends for subsets of diseases, e.g. two IBD phenotypes. We also elucidate a clinically-relevant target COG6 shared between IBD and PSC, as well as a whole range of other target genes missed by the conventional SNP-to-gene assignments methods.

The zebra mussel (Dreissena polymorpha) is an invasive species that causes extensive economic and ecological damage. Here, we identify and characterize the key components of the small RNA (sRNA) and RNA interference (RNAi) pathways in zebra mussels. Like other mollusks, zebra mussels have extensive microRNA (miRNA) and Piwi-interacting RNA (piRNA) machinery but lack or have modified canonical factors needed to produce small interfering RNA (siRNA). Specifically, the zebra mussel Dicer sequence displays substitutions in the conserved DEAD box motif that is required for substrate processivity, and this organism also lacks some attendant accessory factors such as R2D2. We sequenced the small RNA found in both isolated somatic tissue (adductor muscle) and whole animals (including germline), and identified both conserved and novel miRNA and diverse piRNA sequences, but few endogenous siRNAs. To determine whether their remaining sRNA machinery could still be co-opted to initiate gene silencing, we injected dsRNA targeting several genes into zebra mussel adductor muscle. The injected rpn8-targeting dsRNA reduced rpn8 mRNA levels and was processed into sRNA that resemble endogenous miRNAs and piRNAs. The levels of both sRNA types correlated with mRNA knockdown, suggesting that they may act together to initiate RNAi as seen elsewhere. dsRNA targeting other genes produced variable results suggesting that particular criteria may be needed to trigger an RNAi response in this assay. Our results characterize endogenous sRNA pathways in zebra mussels, establish that dsRNA can induce RNAi, and lay the groundwork for further optimizations to establish RNAi-based genetic manipulation tools for this damaging invasive species.

In trypanosmatids genes, transcribed by RNA polymerase II do not have canonical promoters and are organized into directional gene clusters that mature into monocistronic transcripts by a co-transcriptional process known as trans-splicing. Even though gene expression is regulated mainly post-transcriptionally, it is currently understood that chromatin and epigenetics are also involved in this regulation. In eukaryotes, specific signals are normally required for the occurrence of an appropriate transcription initiation. Among them, trimethylation of histone H3 in lysine 4 is the most conserved signal normally detected at transcription start sites of actively transcribed genes. Unlike many model organisms, trypanosomes do not have defined promoters. Instead, transcription initiates in a bidirectional manner from dispersed regions coincident with divergent strand switch regions located between directional gene clusters (DGCs). In T. cruzi, H3K4me3 was observed at the origins of transcription coincident with divergent strand switch regions (dSSRs) in epimastigotes, but it has not been mapped throughout the whole genome at base-pair resolution or in other life stages so far. Here, we set up the CUT&RUN technique for T. cruzi epimastigotes and trypomastigotes. Consistent with a predominant post-transcriptional regulation along the life cycle, we did not find significant differences between life stages. We corroborated that H3K4me3 is enriched at dSSR adjacent to actively expressed DGCs. Moreover, we noticed that this histone mark exhibits different patterns that correlate with the genomic span of the transcription initiation regions and with transcriptional activity. Furthermore, we unveiled that the most actively transcribed DGCs are associated with shorter dSSRs and are located within the core compartment of the genome displaying a more accessible chromatin.

The AMY2B gene encodes pancreatic amylase, a critical enzyme for starch digestion. While previous studies have examined AMY2B copy number variation (CNV) in domestic and some wild animals, less is known about wild carnivores inhabiting regions with limited anthropogenic starch exposure. We analyzed blood samples for serum amylase activity and copy number variation in AMY2B gene from 8 wolves (Canis lupus), 11 brown bears (Ursus arctos), and 3 red foxes (Vulpes vulpes) from Bosnia and Herzegovina. AMY2B gene copy number was assessed using droplet digital PCR (ddPCR), and serum amylase activity and glucose levels were quantified. Although the number of fox samples was limited, foxes and wolves consistently harbored two copies of AMY2B, while brown bears exhibited higher CNV (3.67-8.40, mean 5.88). Serum amylase activity was highest in foxes, moderate in wolves, and variable but lower in bears. Despite differences in AMY2B copy number and serum amylase activity, circulating glucose concentrations did not differ significantly among species. Our findings suggest that variation in AMY2B copy number among wild carnivores may be associated with species-specific evolutionary histories and dietary adaptations, providing insight into genomic mechanisms underlying carbohydrate utilization in natural populations.